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SUMMARY 

The perturbation field induced by a line vortex in a super- 
sonic stream and the downwash behind a supersonic lifting surface 
are examined for the purpose of establishing approximate methods 
for determining the downwash behind supersonic wings . 

Lifting -line methods are presented for calculating super- 
sonic downwash. An unbent lifting line (horseshoe -vortex system) 
is used to compute the downwash behind rectangular and triangular 
wings and the results are compared with the exact linearized solu- 
tions. The chordwise position of the lifting line giving the 
best average agreement with the exact solution is noted for each 
wing. A bent lifting line is used to approximate the triangular 
wing, and the results are in good agreement with the exact solu- 
tion except for points within l/2 chord of the wing trailing edge. 
The use of a bent lifting line seems promising for obtaining 
accurate estimates of the downwash behind swept wings. 


INTRODUCTION 

Several methods for obtaining the downwash behind supersonic 
wings based on linearized theory have been presented. These 
methods utilize conical superposition (reference 1), doublet dis- 
tributions (references 2 and 3), or vortex distributions (refer- 
ences 4 and 5). Each of these methods has certain disadvantages. 
Conical superposition is restricted to wings having plan forms 
composed of straight-line segments and is cumbersome for other 
than trapezoidal or triangular plan foxms. The doublet and vortex 
distributions apply to arbitrary plan forms, but provide integral 
expressions for downwash that are generally very tedious to eval- 
uate. The complexity of these expressions indicates that there is 
a need for a straight-forward procedure for obtaining reasonably 
accurate, if not exact, downwash solutions. 
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A logical approach to the development of a simplified super- 
sonic downwash theory is to derive the supersonic analogues of the 
line -vortex procedures that have proved valuable in subsonic theory . 
Certain differences exist, however, between the properties of 
vortices in a supersonic stream and vortices in a subsonic stream. 
Similarly, the downwash fields behind subsonic and supersonic wings 
differ in certain respects. These differences must be investigated 
before an extension of subsonic techniques is possible . 

The present report, prepared at the NACA Lewis laboratory, has 
three main objectives: (1) The downwash field induced by a super- 

sonic line vortex of constant slope is derived and discussed; 

(2) the downwash behind a supersonic lifting surface is examined 
related to the downwash field induced by line vortices; and 

(3) lifting-line methods for computing downwash are presented and 
calculations based on these methods are compared with the exact 
linearized solutions. Zero -thickness wings (lifting surfaces) are 
considered throughout. 


SYMBOLS 

The following symbols are used throughout this report: 

X,Xi 

y> y i ) Cartesian coordinate system 
z,z t 


u 

V 

w 

i 

T 

k 

I 

i 


perturbation velocities 


unit vectors 


components of vortex vector (three-dimensional field) 
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V 

V 

V 

b 



c 


c 


p 


c e 


. c 


r 


E 


Gi(m) = 


Gi(0) = 

L(y) 

M 


r l = 


U 

Xi = 

Y i = 

Zt = 


components of vortex vector (vortex sheet) 


wing span 

function used in evaluation of finite part of divergent 
integral (equation (7)) 

local wing pressure coefficient 

effective chord of bent lifting line 

root chord of wing 

complete elliptic integral of second kind with modulus 



(Y t - m 1 X 1 )(p 2 m 1 Y 1 - X t ) 
r i [< Y i “ “i x i) 2 + 0- ~ fftnt 2 ^ 2 ] 

-1 *jYj 
r i (Yj 2 + Zj 2 ) 

spanwise lift distribution 
free -stream Mach number 
s] ope of line vortex or lifting line 
fei 2 - 3 2 Yi 2 - p 2 Zi 2 

free -stream velocity (taken in x-direction) 
x -Xi 

y - yi 

i 


z - z 
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angle of attack 

cotangent of Mach angle - 1 

integration interval , 

circulation 

perturbation -velocity potential 

cp T -q> B 

free -stream density 
cross-sectional area of vortex tube 
s panwise distribution of wing circulation 

i 

wing circulation at midspan 

l/2 

resultant vorticity + rf + 

finite part of divergent integral 
line integral along closed curve 
Subscripts j 

a,b points of intersection of forward Mach cone with line vor- 
tex, lifting line, or edge of vortex sheet 

B bottom surface of z - 0 plane 

C line integration 

i = 1,2,3 . . . n points on vortex lines 

l plan-form leading edge 

S surface integration 

T top surface of z = 0 plane 

t plan -form trailing edge 

Y volume integration 

variable of integration 



o 
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Superscripts : 

value of function at point of discontinuity when approach- 
ing from negative y-direction 

value of function at point of discontinuity when approach- 
ing from positive y-direction 


THEORY 

General,, Vortex -Field Relations 

The equations relating velocity and vorticity distributions 
in a slightly perturbed supersonic stream are derived in refer- 
ences 4 and 5. These relations are summarized in the following 
section. The velocity field is assumed to consist of a major 
supersonic free -stream velocity U (taken in the positive 
x-direction) plus small perturbation velocities u, v, a nfl w 
such that the linearized equations of motion are applicable. 

Three-dimensional vortex field. - The velocity field in vector 
form is 



q = (U + u) i + vj + wk (1) 

The vortex vector field (0 is defined as the curl of the velocity 
vector. Thus 


where 


(0 = V X q 

= |i + ti T + £k 


'bw bv' 


( 2 ) 


*) = 


'du dvn 
t3z " 3x J 


. _ /ov _ du\ 
^ \c£c ~ 3y) 
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Vortex lines are lines that are tangent at all points to the local 
vortex vector and are determined from the equation 


dx dy dz 

t = * = t 


(3) 


Reference 4 presents the following integral expressions for 
the perturbation velocities induced by the three-dimensional vortex 
field: 


u 


v 


b 2 I r (y> - T o» 

2 * H V i-o 3 

£ I r (X„S Z„tT 

2 ” J Y 


&x dy dz^ 
o J o o 


dx 0 dy D dz D 


(4a) 


(4b) 


where 


w 



r 

J r 2 

v r ° 


dx 0 dy D dz 0 



(4c) 


r o • V x o 2 - - »v 


The subscript o indicates a variable of integration. The inte- 
gration is conducted over the volume V included in the forward 
Mach cone from the point x,y,z. That is. 




NACA TJJ 1925 


7 


The symbol / designates the finite part of a divergent inte- 

gral. (See references 4, 6, or 7.) The procedure for obtaining 
the finite part is systematized in reference 7 as follows: 

The integrals in equations (4a) to (4c) are of the form 


I 



A(x 0 ) 


(5) 


The upper limit in equation (5) corresponds to limits on the Mach 
cone in equations (4a) to (4c). The finite part of this integral 
is 


where 


(T 



A(x o ) 


(x 2 - x 0 ) 3 /2 




r*i 

= -J(xi) - C 


C = lim 
x o — > x 2 


2A(x 2 ) 
V x 2 “ x o 


J( x o) 


( 6 ) 


(7) 


The term J(x Q ) is the indefinite integral of equation (5) and 
J(xi) is the indefinite integral evaluated at the lower limit. 

The Justification for this procedure is presented in references 4, 
6, and 7. It can be shown that finite parts are treated in a 
manner similar to that for ordinary integrals. The rules of addi- 
tion, differentiation under the integral sign, transformation of 
variables, and integration by parts apply. 

Vortex sheet. - If vorticity exists as a surface of velocity 
discontinuity in the Zq = 0 plane, then outside this plane, |, 

T), and l are all zero, but in the z Q = 0 plane £ = 0 whereas 
| and q are infinite. However, the limits 
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l' - lim | dz 0 

• 4z 0 — >° 


T| ' = lim T] dz c 
dz 0 — >0 

are finite and are given "by (reference 5) 


V = V B " V T 
ti * = u b 


i 


( 8 ) 


The subscripts T and B designate velocities on the top and 
bottom surfaces of the z D = 0 plane, respectively. The perturba- 
tion velocities induced by the vortex sheet are obtained by sub- 
stituting equation (8) in equations (4a) to (4c) . In particular, 
the vertical -perturbation-velocity field (upwash) is given by 


w 



( 9 ) 


The area of integration S includes all the vorticity in the 
forward Mach cone from x,y,z. 


Line vortex. - The vortex lines through all points on an 
.infinitely smal l closed curve bound a vortex "tube." The circu- 
lation 


K = GO O’ 


l/2 

(where ( 1 ) = (£ 2 + t^ 2 + l 2 ) ' is the resultant vorticity and a 
is the infinites imal cross-sectional area of the tube) is constant 
at all points along the tube. The vortex-vector components at 
any point, in tenns of the differential distance dZQ along the 
vortex tube, are 


I 


dx 0 


n 

t 



(*) 


dz 0 

dZ 0 
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The elemental volume is dXodyodZo = odZ Q . A line vortex is gen- 
erated by allowing the cross section of the vortex tube to approach 
zero while maintaining OXJ constant. The upwash induced by such 
a line vortex is found by substituting the preceding expressions 
into equation (4c) and equals 


w = 



( 10 ) 


where the integration is conducted along the portion of the line 
vortex within the forecone from x,y,z. A line vortex cannot 
terminate within a fluid flow field but must either form a closed 
curve or extend to infinity or a boundary of the field. 


Upwash Induced by Line Vortices 

Complicated velocity fields can be generated by the linear 
superposition of relatively simple fields. It will therefore prove 
useful for subsequent developments to determine the upwash field 
induced by line vortices of constant slope. 

Line vortex of constant slope and strength. - The upwash at 
x,y,z due to a line vortex of constant slope m and strength 
K, intersecting the forward Mach cone at Xh.yh (fig. 1) is 
(from equation (10)) 


m 2 (y-mx)dy 0 

[(l-0 2 m?)y o 2 + 2m(0 2 my-x)y o + n^(x 2 -P 2 y 2 -P 2 z 2 )] 3 ^ 2 

(11) 

is performed in appendix A and yields 



The integration 


w = - 


_K_ 

2n 


(y - me) (P 2 my - x) 

,V2 


(x 2 - P 2 ^ 2 - 0 2 z 2 ) 2 [(y - mx) 2 + (1 - 0 2 m 2 )z 2 J 


( 12 ) 
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The finite part of the Integral in equation (11) is obtained by 
substituting the lower limit into the indefinite integral (appen- 
dix A) . No contribution appears from the upper limit . By a trans 
lation of coordinates, the upwash due to a line vortex from xq,yq 
that intersects the forward Mach cone at Xb,7b (fig,. 2(a)) is 
found to be 


where 


w = - 


_K_ 

2rt 


(Yq - mXq)(P 2 mIq - Xq) 


(Yq - mXq) 2 + (1 - P 2 m 2 ) z 2 


(13) 


Xq = X - Xq 

' Yq = 7 - 7q 


r-L =^X X 2 - P 2 Yq 2 - P 2 z 2 

By superposition (fig. 2) of a positive line vortex * from 
xq,yq and a negative vortex -K from y.2>72> the upwash due to 
a line -vortex segment not intersecting the forward Mach cone is 


w = [°2^ ■ M*) 


(14) 


where the notation 


Crq(m) 


*i 


(Yg - m 1 X 1 )(P 2 m 1 Y 1 - Xg) 

(Yg - mqXq) 2 + (1 - P 2 mq 2 ) Zq 2 



(15) 


is used. The subscript for G indicates the appropriate sub- 
scripts for X, Y, Z, r, and m. 

Equation (14) applies for any line vortex of constant strength 
and slope. The circulation K is positive when the vortex vector 
is in the direction of integration. When the line vortex inter- 
sects the forecone, the corresponding limit (infinite Gq(m) term) 
is neglected by application of the finite -part concept. 

Line vortex of constant slope but varying strength. - Line 
vortices may coincide to form a resultant line vortex whose cir- 
culation is the sum of the strengths of the superposed components. 



NACA TN 1925 


11 


The general line vortex is then one of variable strength K along 
the line Xq = x 0 (y 0 ) having the local slope iDq = dy 0 /dx 0 . The 
upwash induced by the segment from Xj,y]_ to ^- 2*^2 W te written 


w 




•c(Y 0 - moX 0 ) 


m 

o o 


dy c 


(16) 


The evaluation of equation (16) is generally tedious. For the 
particular case of constant slope, however, this equation may be 
integrated by parts to yield a useful expression for upwash. Inas- 
much as 

. P 2 < y o ~ PqEq) „ *Qq(°>) 

“c,r 0 3 

the integration by parts gives 



The term K 0 [G 0 (m)] is disregarded at a limit corresponding to a 
point on the forward Mach cone . 


Equation (17) is a generalization of equation (14) to account 
for variations in strength along a line vortex of constant slope. 
Both equations are of fundamental importance because appropriate 
distributions of such line vortices will be used to simulate a 
supersonic lifting surface. 


Characteristics of upwash field due to supersonic line vortex. - 

1. Infinite line vortices of constant strength and slope. An 
infinite line vortex inclined supersonically ( | Pm|-l) to the free 
stream is shown in figure 3(a) . Because both limits in equation (14) 
are neglected, the upwash induced by this line vortex is zero. This 
result agrees with the indications of oblique -airfoil theory because 
the perturbation velocities are zero downstream of the envelope of 
the Mach cones from the trailing edge of a two-dimensional airfoil 
inclined supersonically to the free stream. 
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The Infinite -line vortex along the line (y Q - yq) = m(x Q - xi) 
(fig. 3(h)), inclined euhaanically ( | 0m|<l), has one limit that 
intersects the forward Mach cone whereas the other extends to 
infinity. The upwash is then, for O^0m<l, 


w 


K 

2 « 


lim G 0 (m) 



. (l! - m^) ' 

2* ( Yi . n jx 1 )2 + (l - p 2 m 2 ) z 2 

where is a point on the line vortex. The upwash is infinite 

along the line vortex. For m = 0, equation (18) hecames 


w = 


. — T i . 

2jt Y x 2 + z 2 


(19) 


which is identical with the expression for upwash due to an infinite 
vortex parallel to the stream (along the line y = yq) in an incom- ' 
pressihle field. 


These results indicate that the behavior of the infinite super- 
sonic line vortex for |0m|>l is completely different from that of 
the incompressible flow vortex. However, when |0m|<l, both vortices 
have similar upwash properties in the vicinity of the vortex line 
and are, in fact, identical for m = 0. 

2. Bent line vortices. The upwash due to a bent line vortex 
(fig. 3(c)) of constant strength is 


v - gjj: [Cr^m") - G 1 (m + j] (20) 

where m^** and m^ + designate the slopes of the line vortex before 
and after the bend at xi,y Y . The term r^ appears as a factor in 
the denominator of equation (20) so that the upwash exists only in 
the aftercone from x^yj and is infinite on the cone surface 
(except in the z = 0 plane ) . This infinite value of upwash is 
not to be confused with the infinity introduced at the intersection 
of the line vortex with the forward Mach cone, which is eliminated 
by application of the finite -part concept. 



NACA TN 1925 


13 


Linearized Supersonic -Wing -Theory Relations 

The perturbation velocities on supersonic lifting surfaces 
(zero-thickness wings) have been evaluated for a large variety of 
plan forms . (See, for example, references 7 to 9.) Those results 
will now be utilized to determine the vortex field generated by a 
lifting surface. 

Velocity potential. - If the boundary conditions for a lifting 
surface are specified in the z = 0 plane, the u and v veloc- 
ities are antisymmetric and the w velocities are symmetric about 
this plane. Thus the velocities on the top and bottom surface of 
the z = 0 plane are related by 


( 21 ) 


The discontinuities in u and v constitute a vortex sheet. 
Because the flow is everywhere irrotational, except across this 
sheet, a perturbation-velocity potential cp can be so defined that 

d<p = ^dx + ^dy + ^dz 

= udx + vdy + wdz 

When the undisturbed flow field upstream of the wing is considered 
to be of zero potential, and the boundary condition requiring that 
u ij > = Ujj = 0 off the wing (antisymmetry of u and zero lift off 
wing) is applied, the potential in the z = 0 plane may be . obtained 
by integrating along lines of constant y. 
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where 12 ls the equation of the leading edge as a function of y. 
From equation (23) it may he concluded that 

1. Everywhere in z = 0 plane except behind the wing leading 
edge cp T =q> B = 0. 

2. At a given span station 9^ and 9 b remain constant for 
all values of x downstream of the trailing edge. 

Lines of constant potential for rectangular and triangular wings 
are shown in figure 4. 

Vortex lines. - The equation for the vortex lines (from equa- 
tion (3)) is 

T|* dx - !' dy * 0 (24) 

When the following expressions (equations (8), (21),*and (22)) are 
substituted 

V - ”( V T - v b) 0 - 2 3ST " 2 35T ' 


q 1 = 


o 

U B s 2 5T = 


-2 


ST 


the equation for the vortex lines becomes 


^-dx + ^- dy - 


sr 


dx + 


3F d ^° 


(25) 


Comparison with equation (22) shows that equation (25) represents 
lines of constant potential. Thus the vortex lines coincide with 
the lines of constant potential' in the z = 0 plane . 


Circulation. - The circulation included between two points 
x 1 ,y 1 and X 2 ,y 2 on a wing is given by 


K = 


u dx + v dy + w dz 


(26) 


The path of integration is arbitrary except that the path should 
cross the z = 0 plane only at the two specified points. If the 
integral is taken along the top and bottom surfaces of the z = 0 
plane. 



NACA IN 1925 


15 


52*72 ^i, yi 

K = / (uj dx + vrp dy) + / (ug dx + Vg dy) 

x i>yi ^ 2 ,y 2 


- (^T,2 ” ^T,^ + ^B,l " ^B,2^ = A<P 2 " £t f > l (27) 

where Acp equals <p»p - cpg and represents the jump in potential at 
the point. The quantity Aqp' is, in fact, the doublet potential 
(reference 2), so that the net circulation between two points equals 
the difference in the doublet potential between those points. The 
equivalence of a doublet distribution and a vortex distribution 
indicates that the flow about a lifting surface can be calculated cn 
either basis. 


Circulation aim lift. - The lift per unit span is given by the 
chordwise integration 


L(y) *= ipu 2 


(°p,B * C p,T) dx 


(28) 


After substitution of the linearized values for pressure coefficient 

2u T g 

c p,T B 'T = -jF 

2ub 2 ^ 

°p,b - - — - ' $ 33T 

and integration, equation (28) becomes, because Aq>j = 0, 

L(y) = PU Aq^. (29) 

The factor Acpt is the circulation Included between the leading 
and trailing edges at the s panwise station under consideration. 

When this circulation is designated T, equation (29) becomes 


L(y) = pur 


(30) 
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which is the familiar incompressihle -f low relation. Also 

„r ^ -*B> t 

dy dy 

* (vt “ 

- -d’) t (31) 

Equation (31) relates the shed vorticity to the rate of change of 
spanwise loading. 


APPLICATIONS TO CALCULATIONS OP 
SUPERSONIC DOWNWASH 

The vertical perturbation velocities due to a supersonic line 
vortex of constant slope jhave been presented in a form that permits 
analytical or mschanical evaluation (equations (14) and (17)). The 
vortex distribution associated with a lifting surface has also been 
discussed. These relations will be used to develop exact and approx 
imate methods for calculating downwash behind lifting surfaces at 
supersonic speeds • 


Downwash an Infinite Distance he hind Wing 

The vertical -perturbation -velocity field behind a supersonic 
lifting surface (from equation (9)) is 


w 



(Y o r - X Q V) 


dx 0 


dy Q - 



dx Q dy Q 


(32) 


where the integration over the plan form is designated by S p and 
over the wake by Sy. As x becomes infinite, Xq also becomes 
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infinite in the integral for the hound (plan form) vortices. This 
integral then becomes zero because Xq is of higher order in the 
denominator than in the numerator. Thus, only the Integration over 
the trailing vortex sheet contributes to the vertical perturbation 
velocities at infinity • • The trailing vortex sheet may be con- 
sidered to consist of elemental vortices of infinite length along 
y = constant lines, each having the strength 

4 p 

dK = I' dy Q = - dy 0 . The vertical perturbation velocity at 

®,y,z due to the elemental vortex along y ■ y 0 is, from equa- 
tion (19) 



so that the vertical -perturbation -velocity field at infinity is 
given by 



(33) 


vhere b/2 is the semispan. The velocity field obtained from 
equation (33) is identical to that induced by a subsonic wing with 
the given span loading. The velocity field at infinity is thus 
independent of Mach number (excluding the influence of Mach number 
on dr/dy 0 ) . This result has been derived in references 1 and 2 
by other methods. The evaluation of equation (33) is relatively 
simple and may be used to approximate the downwash several chords 
behind a supersonic wing. 


Regions of Infinite Downwash 

Approximate solutions may modify or introduce singularities 
in the downwash field. It is therefore of interest to establish 
the regions for which linearized theory indicates an infinite 
downwash. 

Infinite downwash in z = 0 plane. - The vertical perturba- 
bation velocities in the z = 0 plane that exist an Infinite 
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distance behind a wing having a discontinuity in dr/dy 0 will 
first he considered. These velocities are determined from the 
equation 


b 



2 


(34) 


The Cauchy principal value is required for points on the vortex 
sheet. For the particular case of triangular loading (fig. 5(a)) 


dr 2l m 

dy Q ” t 


in the integration interval 


T- y o-° and 


dr _ zr m 

4y 0 “ ” * 


in the interval 



(r m is the circulation at the midspan) . 


When these values for 
integral yields 


dT 


are substituted in equation (34), the 



(35) 


Infinite upwash exists along the lines y = ±b/2 whereas infinite 
downvash exists along the line y = 0. These infinite values are 
due to the discontinuities in the spanwlse vorticity distribution 
nr>a apply for all points on these lines downstream of the wing 
trailing edge. 

In general, if (dr/dy 0 )“ and (dT/dy 0 ) + represent a dis- 
continuity in the rate of change of spanwlse loading at station 
y = y^, infinite vertical perturbation velocities, will exist along 
y = y, downstream of the trailing edge. For (dr/dy 0 ) “ < (dT/dy 0 ) + 
infinite upwash will exist, and for (dr/dy 0 )“ >(dr/dy 0 ) + infinite 
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dovnwash will exist along this line. Such discontinuities in 
ar/dyo originate both at a wing tip and at those points along a 
supersonic trailing edge where the plan-form slope is discontinuous 
and the local -wing -pressure coefficient is not equal to zero. This 
discontinuity in vorticity may he verified by the methods employed 
in reference 1 for finding the upwash and sldewash directly behind 
a supersonic trailing edge. 

The discontinuity in shed vorticity at the tips of an ellip- 
tically loaded wing is a special case of the previously stated rule 
governing infinite vertical perturbation velocities in the z = 0 
plane. For wing loading given by 

r = “ ? 7 ° 2 


the shed vorticity is 


dF 

dTo = 



^m 7o 




and the vertical -perturbation-velocity field at infinity, in the 
z = 0 plane, is 



Integration yields 



for 
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and. 



As indicated in figure 5(b), the vertical perturbation velocity is 
discontinuous at y ■ ±b/2, but is bounded for all points on the 
vortex sheet . 

Infinite downwash on Mach cone from wing tip. - Reference 3 
indicates infinite vertical perturbation velocities on the down- 
stream Mach cones from the tips of a triangular wing. This result 
will now be extended to apply to any wing tip formed by the inter- 
section of a subsonic leading edge and a supersonic trailing edge 
provided that the slope of the subsonic edge is not zero at the 
tip (fig. 6). 

The contribution to the vertical -perturbation -velocity field 
due to the bending of an elemental vortex at the trailing edge 
(fig. 6) is, from equation (20), 

— dy 

dw - - Ay ° 2 - - ° [o 0 (m")- - G o (0)] (36) 


where 


G o (0) = ± 


Vo 


•o (Y c 


+ Zf 


The vertical perturbation velocity at a point on the Mach cone from 
this tip due to the bending of the vortices is found by integrating 
(along the trailing edge) 


w 



G 0 (m ) " G o (0) 



d y 0 


(37) 
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Equation (37) in the expanded form becomes 

-if 

J 7a L 


(y 0 - yx 0 )(p 2 yr 0 - x 0 ) w _ w 

(Y 0 - V x o) 2 + (1 - p 2 m 0 " 2 )z 2 + Y q 2 + z 2 


*o*o 


i dr . 

X TT f 


* 0*^0 


(38) 


The limits of Integration are roots of r Q , so that the factor 
l/r 0 is singular of order l/2 at the limits and the integral is 
improper (assuming (m^, )y^ . g / 0) , The convergence of the Integral 

depends on the nature of dr/dyo at the limits. Appendix B shows, 
however, that dT/dy 0 is also singular of order l/2 at the tip. 

The combined singularity causes equation (38) to diverge at the 
upper limit and results in logarithmically infinite vertical per- 
turbation velocities on the Mach cone from the tip. The divergence 
is a consequence of both the singularity in dl/dy 0 and the singu- 
larity on the Mach cone from a bend in the elemental line vortex. 


The infinite vertical perturbation velocities on the Mach cone 
from the tip do not appear in the z = 0 plane inasmuch as equa- 
tion (38) then reduces to 


v 


i_ f V ko - fjV ar , 

2* J_ (T 0 - Bo-Xo)T 0 Uo ° 


and the singularity due to r Q is no longer present. 


Approximate Downwash Solutions 

Several approximate methods for obtaining downwash were con- 
sidered. Methods based on a lifting line seem the most promising 
and are discussed in the following sections. A bent lifting line 
is proposed for determining the downwash behind a swept wing and 
an unbent lifting line (horseshoe vortex system) for determining 
the downwash behind an unswept wing. These methods are applied to 
compute the downwash behind triangular and rectangular wings 
the results are compared with the exact linearized solutions . 
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A lifting line concentrates the chordwise loading into a line. 
Thus the hound circulation is represented by a line vortex of vari- 
able strength (k = T = &p>t), whereas the t railing-vortex system 
maintains the sane vorticity (£' - -dT/dy 0 ) hut now originates at 
the line rather than at the trailing edge. 


Bent lifting line. - A lifting line approximating the section 
centers of pressure 


<Cp,B - C p> t)(x - xj) dx 


‘op - + 


/ 


(°P,B “ C p,T^ ** 


se ems to he a reasonable representation of a sweptback or swept- 
forward wing. In order to facilitate downwash calculations, the 
line of section centers of pressure can he approximated hy two 
straight-line segments, each connecting an end point to the mid- 
point of the line. The result is a bent lifting line (fig. 7) of 
span h and effective chord c e (x-distance between midpoint and 
end points of lifting line) . For a sweptback lifting line, the 
slopes at the midpoint are 

( m o")y 0 = 0 = " 2c^ 

and 

Wy 0 - ° - zt; 

The vertical perturbation velocities induced by the bound 
vortices can be determined from equation (17) and equal 


G ° (m) ty^ (39) 

where the integration is conducted along the lifting line (figs. 8(a) 
and 8(b)). The vertical perturbation velocities induced by the 
trailing -vortex system are obtained from 
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The term G o (0) is as defined for equation (36) and the integration 
is again conducted along the lifting line. The vertical perturba- 
tion velocities induced by the complete lifting-line system (bound 
and trailing vortices) are then the sum of equations (39) and (40) 
and equal 




The value = (m 0 ")y 0 = o 


is used for the integration interval 


7a— 7q — Q ®- nd m o = (m r + )_ _ n is used for the interval 

o«y 0 Sy b . ,0 " 


Equation (41) is the expression for the vertical perturbation 
velocities behind a swept wing using the bent -lifting -line approxi- 
mation. This equation can also be considered as derived from a 
superposition of a system of bent line vortices of constant strength 
(equation (20) ) , as indicated in figure 9. The nonintegral term of 
equation (41) is then the contribution of the bent line vortex of 
strength r m whereas the integral term represents the contribution 

of the elemental line vortices of strength “ dy rt . 

d y 0 


The, integral term of equation (41) may be evaluated by analyt- 
ical or mechanical methods. When mechanical methods are used, the 
singularities in the integrand must be isolated. Suitable proce- 
dures for isolating the commonly encountered singularities are as 
follows : 


(1) Singularity due to intersection of forward Mach cone with 
lifting line. The integrand in equation (41) is infinite at the 
intersection of the forward Mach cone and the lifting line (for 
z / 0), The contribution 5w of the integral, for the interval 
yi) - eb<y 0 <yi 3 (fig. 8(a)), to the vertical perturbation velocity 
field may be written 





\ 
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for 

Pm Q >l, 



where m Q is the slope of the lifting line at y^. A s imilar pro- 
cedure applies for the singularity at y a . 

(2) Singularity in dr/dy c at wing tip. The vertical pertur- 
bation velocity at a point whose forward Mach cone intersects the 
edges of the trailing vortex sheet (fig. 8(b)) is obtained from '• 

v v 

equation (41) by integration between the limits - g<y 0 <g. If 

dT/dy 0 is singular at these limits, the singularity may be isolated 
by a procedure similar to that used in equation (42). Thus, the 
contribution to the vertical perturbation velocity field, from the 

■u n. 

integration interval g- - < y Q < g-, may be written 
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(5) Singularity at points on vortex sheet. The Cauchy prin- 
cipal value of the integral (equation (41)) is required for deter- 
mining the vertical perturbation velocity at a point on the vortex 

sheet z = 0, - |^y<|. For this case, considering the interval 
y - €<y <y + e (fig* 8(c)), the contribution of the integral to 

the vertical -perturbation-velocity field is 


Sw 



■» 0 ^0 2 - eV ar 
(y c - *oX 0 )y 0 <*y 0 ° 


( 44 ) 


If dr/dy© can be approximated by the first two terms of the 
Taylor's expansion 
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= y 


+ (y 0 - y) 



and c is sufficiently small that 


■W x o 2 - - - 1 

(Y 0 - 

equation (44) may "be written 



Inasmuch as &T/dy 0 is an odd function in y 0 for a wing symmet- 
rical about the midspan, dr/dy Q is either discontinuous or zero 
for y D = 0. If dr/dy 0 is discontinuous, the vertical perturba- 
tion velocities are infinite along this line. 

Unbent lifting line. - The unbent lifting line ( horseshoe - 
vortex system) appears to be a reasonable representation for an 
upswept wing. The use of such a lifting line considerably sim- 
plifies the calculations. 


The vertical perturbation velocities induced by an unbent 
lifting line (along the y-axis) are obtained from equation (41) 
by setting eiq = ® and are equal to 


w 



x Y 0 (r c 2 - p 2 z 2 ) aT_ 
r Q (x 2 - 0 2 z 2 )(Y q 2 + z 2 ) d-y o y ° 


( 46 ) 


The singularities are isolated as follows : 



28 


NACA TN 1925 


(1) Singularity due to intersection of forward Mach cone with 
lifting line . 


5w = -sr 





- 


f 

x Y 0 (r 0 2 - P 2 z 2 ) 

dr 

( p 2 z 2 dr \ 

d y Q 

j 

(x 2 - p 2 z 2 )(Y 0 2 + z 2 ) 

*y 0 + 

\x y 0 d y 0 / 

r o 


_ 


y o=yi> 



A. ( ££ 4L.' 

2 « \x Yq d y 0 , 


yo^h 






d y c 


where 


7b 


d y 0 2 , -l / e b 
- = ^ 8in 


*o P 


'y-b _e t 


(2) Singularity in df/dy 0 at the tip. 


*T> 

2"S 


x Y Q (r o 2 -p 2 z 2 ) 


r Q (x 2 -P 2 z 2 )(Y 0 2 +z 2 ) 


x Y 0 (r 0 2 -P 2 z 2 ) 


r 0 (z 2 -P 2 z 2 )(Y 0 2 +z 2 ) 




J h 
T o=2 


1_ 

2rt 


x Y 0 (r o 2 - P 2 z 2 ) 


rjx 2 - P 2 z 2 )(Y q 2 + z 2 ) 


y o -2 


(r) b r 

b y o = 2 - C b 


(3)- Singularity at points on vortex sheet. 

^d 2 r 


5w = l| 


y 0 - y 
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The ch'ordwlse location of the unbent lifting line that will give 
the best average agreement with the exact linearized solution is 
still to be determined. According to the techniques used in subsonic 
wing theory, an unbent lifting line at the wing center of pressure 
should be a good first approximation. However, further investi- 
gation is required. It may be possible to determine the best loca- 
tion for each general class of plan forms by comparing the lifting- 
line solution with the exact linearized solution for several 
representative plan forms. 

Examples and discussion. - A bent lifting line and an unbent 
lifting line will be used to obtain solutions for the downvash in 
the y = 0 plane behind triangular wings having subsonic leading 
edges. The chordwise distribution of wing loading suggests the use 
of the bent lifting line . The unbent lifting line will be used for 
purposes of comparison. An unbent lifting line also will be used 
to compute the downwash in the y = 0 plane behind rectangular 
wings. The chordwise location of the unbent lifting line giving 
the best average agreement with the exact linearized solutions 
will be noted for both the triangular and the rectangular wings. 

In order to simplify the expressions, M = ,./2 (that is. 6 = 1) 
will be assumed. ' ' 


The spanwise circulation distribution for a triangular wing 

of span b and root chord c r , having subsonic (~-<lj leading 
edges, is (from equation (24) of reference 2) ' ° r ' 

T = ^ /T TXT 2 

e y- ^2 7o 

where E is t he complete elliptic integral of the second kind 


with modulus 


V 1 - (£)’ • 


The circulation at the mids pan is 
P _ cffjb 

IT 


and the rate of change of circulation is 


dr 4au y G 
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The expression for the vertical perturbation velocities in the 
y = 0 plane using the bent -lifting-line approximation of fig- 
ure 10(a) is 


w 



[G 0 (m+) - G- o (0)] 
. /l - 4r y. 2 


Vn <&r 


(47) 


vhere (m 0 + )^ = Q 


— . Equation (47) is obtained from equation (41) 

c r 


ar>a the properties of an even function. The vertical -perturbation- 
velocity field behind an unbent lifting line , having the same 
loading (fig. 10(b)) is from equation (46), 



The 


ratio 


a panwise circulation for a rectangular wing of aspect 

^ Td 

for the tip region ^ - c r ^ y 0 is 


i>2 



and for 



(49a) 


T = 2cflJ c (49b) 

r 

Equations (49a) and (49b) were obtained from equation (20) of refer- 
ence 8. The wing is illustrated in figure 10(c). The rate of 

change of spanwise circulation for - c r < y Q < ^ 
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ar 



The vertical -perturbation-velocity field in the y = 0 plane 
behind an unbent lifting line having this loading is 


w 



x (r Q 2 - z 2 ) y Q 

- z 2 )(y o 2 + z 2 ) 



d 7o 


(50) 


Equation (47) was evaluated to determine the downwash along 
the line y = 0, z = 0 and along the line y = 0, z = 0.1 b 
for triangular wings with aspect ratios (2b/c r ) of 1.6 and 3.2. 
Equation (48) was evaluated for the same wings at points along 
the line y = 0, z = 0. Equation (50) was evaluated at points 
along the line y = 0, z = 0 for rectangular wings with aspect 
ratios (b/c r ) of 2 and 4. The integrations were performed mechan- 
ically for equations (47) and (50) and analytically for equa- 
tion (48) . The results are compared with the exact linearized 
solutions obtained from references 1 to 3 in figures 11 to 14. 


The bent -lifting-line solutions for the triangular wings are 
shown in figures 11 and 12. The discontinuity in the curves for 
the z = 0.1 b plane (figs. 11(b) and 12(b)) designates large 
negative (upwash) values that become infinite on the Mach cone 
from the tip. as indicated in the discussion associated with 
equation (38). The agreement with the exact solutions is good 
for all points except those within l/2 chord from the trailing 
edge. The correlation is unexpected because of the large contri- 
bution of the nonintegral term of equation (47). This term is 
associated with the bending of the lifting line at the midpoint. 
This bend was artificially introduced. The agreement with the 
exact downwash solutions, however, indicates that the bent lifting 
line is a good average representation of the triangular-type wing 
(at least in regard to the downwash solution in the region of 
the line y = 0, z = 0) . This method should give even better 
agreement when used to represent a sweptbac!'. wing such as the one 
represented in figure 7, because the bent lifting line would then 
more closely approximate the actual vorticity distribution. 
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The uribent-lifting-line solution for the triangular wings is 
presented in figure 13. The lifting line is placed at the 
3/4-chord point to give the test average agreement with the exact 
linearized solution. (The center of pressure for these wings is 
at the 2/3-chord point.) The agreement is not as good as that 
obtained with the bent lifting line and indicates that the bent 
line is more suitable for computing downwash behind triangular 
wings. The smaller-aspect -ratio wing (2b/c r = 1.6) is in better 
agreement with the exact linearized solution because of the rapidity 
with which the downwash approaches the asymptotic value at infinity . 

The unbent -lifting-line solution for the rectangular wings is 
presented in figure 14. The downwash obtained with the lifting 
line at the l/2 -chord point is in excellent agreement with the 
exact linearized solution. (The location of the. center of pressure 
is at the 4/9-chord and 10/21-chord points for aspect ratios of 2 
and 4, respectively.) These calculations suggest that the unbent 
lifting line is well suited for computing the downwash behind 
unswept supersonic wings and that the best chordwise position for 
the unbent lifting line is at or slightly downstream of the wing 
center of pressure. 

A more accurate estimate of the downwash in the z = 0 plane 
close to the trailing edge may be obtained by Judiciously fairing 
the curve obtained by the line -vortex method to the known value 
of downwash at the trailing edge. At a subsonic trailing edge 
satisfying the Kutta condition, -w/aU is unity; whereas at a 
supersonic trailing edge, -w/cUJ can be computed by the method 
presented in reference 1. 

It should be noted that linearized theory neglects the effect 
on the downwash of the friction wake and the displacement and dis- 
tortion of the trailing vortex sheet. Aa experimental program, 
such as that reported in reference 10 for wings in subsonic flight, 
is ultimately required to determine the necessary modifications to 
linearized theory that will result in good agreement between theory 
and practice. 


. SUMMARY OF ANALYSIS AND APPLICATIONS 

The perturbation field induced by a line vortex in a super- 
sonic stre am and the downwash field behind a supersonic lifting 
surface have been examined for the purpose of establishing approx 
j mat e methods for the calculation of supersonic downwash. 
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An Infinite line vortex of constant strength and slope, super- 
sonically inclined to the free stream, induces no perturbation 
field. A suhsonically inclined line vortex has properties similar 
to those of a vortex in an incompressible -flow field. Bends in a 
line vortex induce infinite vertical perturbation velocities on 
the surface of the downstream Mach cone from the bend (except in 
the z = 0 plane ) . 


The downwash field behind a supersonic lifting surface differs 
from that behind a subsonic wing in several respects. For a super- 
sonic lifting surface, discontinuities in shed vorticity occur at 
those points along a supersonic trailing edge where the plan-form 
slope is discontinuous and the local pressure coefficient is not 
zero. These discontinuities lead to singularities in the downwash 
field in the z = 0 plane. Also, the vertical perturbation veloc- 
ities are logarithmically infinite on the downstream Mach cone 
from a wing tip formed by the intersection of a subsonic leading 
plan -form edge and a supersonic trailing plan-form edge. 


A bent -lifting -line method has been proposed for the solution 


ox one uuwiiWaSh field behind swept wings . 


When applied to a 


triangular wing, this method gave results that were in very good 
agreement with the exact linearized solution for points near the 
line y => 0, z = 0 except for points within l/2 chord of the 
wing trailing edge. 


An unbent lifting line (horse shoe -vortex system) haa been 
proposed for unswept wings. This method was applied to determine 
the downwash behind rectangular wings with aspect ratios b/c r of 
2 and 4. Excellent agreement was obtained for both aspect ratios 
by placing the lifting line at the l/2 -chord point. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, April 29, 1949. 
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APPENDIX A 


UPWASH INDUCED BY UNE VORTEX 

The upwash Induced "by the line vortex segment of figure 1 is, 


3 2 K 
" 2n 


I C • m 2 (y-mx)dy 0 

J 0 (i-p 2 n£)y 0 2 + 2m(P 2 my-x)y 0 + n^(x 2 -P 2 y 2 -3 2 z 2 ) ^ 


(ID 

The upper limit is at the intersection of the line vortex (y 0 = mx D ) 
with the trace of the forward Mach cone in the z = 0 plane 


k- 


x )2 _ p 2 (y - y 0 )' <i - p 2 z 2 = oj and is therefore the appropriate 
root of°the expression appearing in the denominator of the integrand. 
The roots are 




-m( E 2 mv - x) ± Pm ^y - mx) 2 + (1 - P 2 m 2 )z^ 
y a^h = (i I ($2 m 2) 


so that equation (11) may "be rewritten 


(Al) 



-2it(l - P 2 m 2 ) 2 |r II 

P 2 m 2 K(y - mx) ~ U [fra “ Jo)^ " ^o)] 3 ' 2 

From equations (6) and (7), 

• (T = - j(o) - c 


(A2) 


(A3) 


where 


2(y a + Tb) 
(y a - y h) 2 


J(0) 
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and 


C = lim 
?o — ^b 


= lim 

y 0 — > y-b 

= 0 


2 A (yb) T , v 

V^b - y 0 7 ° 


2(y a + ?b - 2y 0 ) 


[fra-y^ 


(y a -y t ) 2 y^a-y^Hyb-y^ 


Thus 


/i- - 


2(y a + yb) 


(A4) 


Substituting equation (A4) in equation (A2) and solving for w 
yields 


v = - 


_K_ 

2rt 


(y-ng)(P 2 my-x) 


J. 

/ 2 n 2 2 a 2 2 v 2 |\ x2 ,, „2 2, 2 

[x -P y -0 z ) [(y-mx) + (1-P m )z 


(A5) 
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APPENDIX B 


LOADING IN VICINITY OF WING TIP 

The nature of the loading in the vicinity of a wing tip formed 
by the intersection of a subsonic leading edge with a supersonic 
trailing edge (fig. 6) is to be determined. 

According to linearized theory, uj in the region of a sub- 
sonic leading edge is singular of the form 


«T 



where xj is the equation for the leading edge as a function of 
y Q . This relation may be deduced from equation (11) of refer- 
ence 10. The corresponding wing circulation is 


I dx Q 

r = Acp + o 2 / — = 4 ^/x^. - Xj (B1 ) 

/ V*o - *7 

x i 


The derivative of equation (Bl) is 


dr 

^o 


d(Xfc - x z ) 
dy 0 
V* t - 


(B2) 


where (xt - xj) is the wing chord as a function of y 0 . However, 
(xt - must be of the form 


*t ‘ X I • y 0 C f(y 0 S (B3) 


(where jf(y Q )] t Q»°°) in order to satisfy the restrictions 

7 0 “ 0 

that at the wing tip 
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(a) the chord is zero 


(X * ' = o - 0 

(b) the slopes of the leading and trailing edges are neither 
equal 

d <*t “ X Z } 

d y 0 _ 

w 17 u 

nor in the free-stream direction 



~&i*t ~ x z) 

dy o . 


Substituting equation (B3) into the denominator of equation (B2) 
yields 



2 d(x t ~ X Z ) 
ar = d y 0 
dy o [ f (y 0 )] 


(B4) 


Equation (B4) indicates that dr/dy Q is singular of order l/z at 
y Q = 0 for the wing tip of figure 6. 
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Figure 1. - Geometric relations for determination of upwash 

Induced by line vortex. 
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Line vortex of 
strength K 


W 1 = ' A G l (m) 


(a) Line vortex from (x 1# y 1 ). 




Figure 2* - Superposition for obtaining upwash induced by line- 
vortex segment from x^ # y]. to 
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(x,y,z) 


(a) Infinite line vortex inclined supersonically 
to free stream. 



(x,y,z) 


fb) Infinite line vortex inclined subsonically 
to free stream* 



Figure 3* - Typical line vortices 



NACA TN 1925. 



(a) Rectangular wing. 



Figure 4. - Lines of constant potential for rectangular and 

triangular wings. 
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(a) Intersection of forecone with lifting line. 



(b) Intersection of forecone with edge of vortex sheet. 



Figure 8. - Inproper intervals of dcwnwash integral. 
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(a) Bent-lifting-line representation o<f 
triangular wing. 



(b) Uribent-lif ting-line representation 
of triangular wing. 



Figure 10. - Lifting-line representations of triangular and rectangular 
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Exact linearized solution 
(reference 2) 

Bent-lifting-line solution 
(equation (47)) 
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Figure 12. - Downwash in y= 0 plane behind triangular wing with aspect ratio of 3.2 

using bent lifting line. M = Vg*. c r 
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